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Abstract 

We derive the angular power spectra of intensity and polarization of Galactic syn- 
chrotron emission in the range 36 < £ ^ 10^ from the Parkes survey mapping the 
southern Galactic plane at 2.4 GHz. The polarization spectra of both electric and 
magnetic parity up to £ ~ 10^ are approximated very well by power laws with slope 
coefficients ~ 1.4, quite different from the CMB spectra. We show that no problem 
should arise from Galactic synchrotron for measurements of CMB polarization in 
the cosmological window. 
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1 Introduction and main results 

The purpose of this work is to compute the angular power spectra of syn- 
chrotron emission from the Galaxy, and in particular of its polarized com- 
ponent. Synchrotron radiation is expected to dominate the linearly polarized 
component of the sky background in a wide frequency range up to a few 
tens of GHz. It is important for the knowledge of the Galactic structure, and 
also because it will be a major contaminant in forthcoming measurements of 
the polarization of the cosmic microwave background (CMB), planned from 
ground (Keating et al., 1998; Sironi et al., 1997; Pisano, 1999; Hedman, 1999) 
and from space (Wright, 1999; Cortiglioni et al., 1999; De Zotti et al., 1999). 

Separation of the various contributions to the total background in the "cos- 
mological window" (say, ~ 50 — 90 GHz) promises to be a dehcate job, and 
its success will rely on sensible assumptions about the spectral and spatial 
properties of Galactic foregrounds (e.g., Tegmark et al. 2000, Prunet et al. 
2000). In this connection we observe that at frequencies of tens of GHz the 
dominant contribution to the angular spectrum may be a Galactic foreground 
(synchrotron or dust), the integrated background of unresolved extragalactic 
sources, the primordial anisotropy of CMB or a non-primordial anisotropy (for 
instance arising from the gravitational lens effect) depending on the range of 
the spherical-harmonic index /. Therefore the cosmological window should be 
understood as a bidimensional region of the {u, I) plane, which is furthermore 
different for anisotropy and polarization. For this reason angular power spec- 
tra, which provide a standard tool for CMB analysis, are becoming relatively 
commonplace also for foregrounds, although in the latter case they provide 
an incomplete description because of the strong phase coherence of Galactic 
emission. Intensity (or temperature) power spectra of foregrounds are reason- 
ably approximated by power laws, Cn oc with slopes steeper than CMB. 
For the latter the scale-invariant adiabatic-wave spectrum predicts af^^ ~ 2 
on large scales and power excess at 10^ ^ £ ^ 10^ because of several acoustic 
peaks. Both the large-scale behaviour and the first acoustic peak appear to be 
confirmed by experiment (Barreiro, 1999, de Bernardis et al. 2000, Hanany et 
al. 2000). For synchrotron, values af^ ~ 2.5 3 have been derived from the 
408-MHz maps of Haslam et al. (1981) by Tegmark & Efstathiou (1996) and 
Bouchet et al. (1996), and a^^^ — 3 from the 1420-MHz northern sky survey 
(Reich & Reich, 1986) by Bouchet & Gispert (1999). However, the analysis of 
Tenerife patch at the same frequencies provides a nearly scale-invariant spec- 
trum, aY"^ ~ 2, except near the resolution limit of the maps (Lasenby, 1997). 
For dust emission Gautier et al. (1992) gave aj"^^ ~ 3 down to the IRAS 
resolution of a few arc minutes, and a similar result was found for DIRBE 
(Wright, 1998); from the combined DIRBE and IRAS maps Schlegel et al. 
(1998) derive af'^^^ ~ 2.5. The situation is less clear for free-free emission: 
While Kogut et al. (1996) find a^^ ~ 3 correlating COBE-DMR with DIRBE, 
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Veeraraghavan & Davies (1997) derive aj^ ^ 2.3 at 53 GHz from Ha maps but 
with a much lower normalization. The discrepancy supports the case for dust 
emission in the cosmological window, the mechanism of rotational excitation 
of small grains having been proposed by Draine & Lazarian (1998). 

Angular spectra of polarized emission have been measured so far neither 
for CMB nor for Galactic foregrounds. Theoretical modeling for CMB gives 
(for the electric-parity polarization excited by density waves in the standard- 
model) a^'^^ rather close to zero, which means much smaller spatial correla- 
tions at large angular scales than for temperature fluctuations. A correspond- 
ing excess of small-scale structure (in comparison to temperature fluctuations) 
is found by Prunet et al. (1998) for polarized dust emission in modeling based 
on the HI maps of the Leiden-Dwingloo survey [see also Sethi et al. (1998) 
and Prunet & Lazarian (1999)]: Their results are a^^^ ~ 1.3 and aj^^^ ~ 1.4 
for the electric and magnetic parity spectra, respectively. The angular power 
spectrum is therefore steeper for Galactic dust than for CMB. 

Now an important question should be posed: can this result be generalized 
to other foregrounds? If so, the separation of CMB polarization signals would 
be easier. The present work considers the polarization spectra of Galactic 
synchrotron, and can be regarded as complementary to that of Prunet et 
al. (1998) since polarized free-free emission should not be dominant at any 
frequency. Our approach is the analysis of low-frequency maps. Unfortunately 
no full-sky survey is available for the polarized component of synchrotron, the 
largest coverage being provided by Brouw & Spoelstra (1976) at the expense 
of a quite sparse sampling. Here we make use of the Parkes survey of the 
southern Galactic plane (Duncan et al., 1995, 1997) which, although limited 
to a relatively small window of ~ 0.39 sr, provides uniformly covered maps, 
not suffering from undersampling. Its frequency of 2.4 GHz probably allows 
a reasonable extrapolation to the cosmological window, thanks to the limited 
effects expected from Faraday rotation. The typical values of 5-10 rad m~^ 
for the rotation measure (RM) reported by Spoelstra (1984) and Wieringa et 
al (1993) allow to estimate a negligible depolarization of few percent at 2.4 
GHz. 

The main result of this paper is that the synchrotron polarization spectra up 
to £ ~ 10^ are approximated very well by power laws with slope coefficients 
oe — ctB — 1-4, a result very close to that of Prunet et al. (1998). The intensity 
spectrum is also found to have a moderate slope, close to that of polarization 
spectra, unless we use alternative maps (available at the Parkes WEB siteQ) 
where Galactic-plane sources have been removed; therefore this result should 
not be extrapolated to the whole sky. In the last section of this paper we 
provide reasons, suggesting that the slopes of polarization spectra are likely to 
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apply to the whole sky. We show that if this assumption is correct, synchrotron 
emission cannot be a problem for CMB polarization measurements in the 
cosmological window. Otherwise, a value ~ 1.4 can be used as a lower limit to 
ue and a^; then the slope differences with respect to the CMB polarization 
spectra would still be larger, and the CMB polarized signal would prevail even 
more strongly in the cosmological window. 



2 Angulair spectra in the Parkes survey 



Figure 1 shows the sky coverage of the Parkes 2.417 GHz survey and other 
available surveys as well. The Parkes survey covers a strip, 127° long and at 
least 10° wide, with a FWHM resolution of 10'. 4. Maps are available at the 
Parkes WEB site for the total intensity / both before and after subtraction 
of structures with size < 2°, and for the Stokes parameters Q and U with no 
source subtraction. 

Because of the limited sky coverage, spherical-harmonic spectra are suitably 
obtained from a standard Fourier analysis according to the technique devel- 
oped by Seljak (1997) for the CMB polarization and applied by Prunet et 
al. (1998) to the polarized dust foreground. After the Fourier components of 
Stokes parameters /(I). Q{1) and U(l) are computed on a map covering a 
solid angle Q, the estimators for the power spectra of intensity Cj£ and total 
polarization Cp^ can be readily computed by means of the equations 



E - ^7'} b-'{i), (1) 

cp,^\^§-^j:[Qm%i)+um*{i)]-2w-p'Y-\e), (2) 

where the sums are performed over the modes with wavevector magni- 
tude around i, b{£) is the window function and wjp = ^laj p/A^pixci is the 
pixel-independent measure of noise with aj^p the pixel noise. The electric and 
magnetic parity of polarization can also be separately computed by means of 

^i^^ = 1^ E 1^(1) cos(20i) + ^(1) sin(20,)|^ - ^pij 6-2(£), (3) 
= E l-Q(l) «in(20i) + U{\) cos(2</.,)|^ - ^pij b-\i), (4) 
with 01 the direction angle of 1. 
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Fig. 1. Sky coverages of synchrotron polarization surveys. We report the regions 
observed at 0.408-1.411 GHz by Brouw & Spoelstra (1976); in the 2.417-GHz Parkes 
survey (Duncan et al., 1995, 1997); at Effelsberg, 2.695 GHz (Junkes et al. 1987; 
Duncan et al., 1999) and 1.4 GHz (Uyaniker et al. 1999); and in the Canadian 
Galactic Plane survey (CGPS, 0.408 and 1.420 GHz) (English et al., 1998). 

For the Parkes survey the window function can be approximated by a Gaus- 
sian, b{£) = exp + l)(T^/2] with cxb = 4'. 4 . The polarization noise is not 
constant, being 11 mJy beam area~^ for most of the sky coverage but as low 
as 6 mJy beam area"-*^ in some regions. We chose to extract from the survey 
constant-noise square patches of 10° x 10° (i.e., 150 x 150 pixels). Wc thereby 
obtained 12 independent submaps (the first one being centered at / = 360° 
and moving toward decreasing Galactic longitudes), and on each of them we 
performed the Fourier analysis and derived a set of four spectra Cxe (with 
X = I, P, E and B) by means of Eqs. (l)-(4) (the P spectrum was in fact re- 
dundant, but we computed it for all of the submaps for checks of consistency). 
A potential problem with the Fourier analysis, when it is applied to small 
maps, arises from border effects, which may give spurious contributions to 
the spectra. A solution to this problem is considered by Hobson & Magueijo 
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Fig. 2. The input and output CMB spectra (dotted and full line, respectively) for 
a test of Fourier analysis in the Parkes survey window. The input spectrum is 
computed for a standard CDM model with a secondary-ionization optical depth 

Tion = 0.18. 

(1996), who suggest to modulate the signal by some function vanishing at 
the borders (apodization). We carried out our analysis both with appropri- 
ate cosine functions and with no modulation. The differences turned out to 
be small in the range i ^ 1000, except for the intensity maps with source 
subtraction. In Figures 3-5 we explicitly report the results obtained with the 
cosine modulation up to £ = 1100, which is close to the resolution limit of the 
survey. In order to test the reliability of the computed spectra, we also built 
up simulated maps for CMB, using angular spectra computed by means of 
CMBFAST (Seljak & Zaldarriaga, 1996) and HEALPix package^, and per- 
formed the Fourier analysis on them within the Parkes survey window. The 
results shown in Figure 2 prove that the method is quite reliable: The output 
spectrum differs from the input one only for some smoothing of the sharpest 
features, which is particularly evident for low values of i. Smoother spectra - 
such as those expected for synchrotron - are less affected by the procedure. 
The results for all of the 12 submaps also are mutually consistent in the CMB 
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Fig. 3. The angular power spectra Cn and Ce(. (dashed and full lines, respectively) 
for the 12 submaps of the Parkes survey. 

case. 

The spectra obtained from synchrotron submaps, however, show large varia- 
tions in the normahzation; this fact is certainly not unexpected, due to the 
phase coherence of Galactic structure. Figure 3 shows the results for Cj^ (with 
no source subtraction) and Cei. from the individual submaps; the results for 
magnetic parity polarization are similar to those for electric parity (the in- 
tensity spectra after source removal, however, are much steeper than those 
reported in the Figure, as we will discuss below). 

The normalization variations for intensity and polarization cover different 
ranges (spanning more than 3 and about 2 orders of magnitude, respectively), 
and they are not strictly correlated since the polarization degree is far from be- 
ing constant. The slopes of the curves up to £ 10^, however, are much more 
mutually consistent than normalization. Visual inspection of Fig. 3 shows a 
couple of low-emission regions (submaps 9 and 12, placed at / < 280°) with 
particularly flat intensity curves at £ ^ 10^, but this feature does not appear 
in polarized emission. Clearly most of the curves are reasonably approximated 
by power laws, but when we fit the data with functions Cxt — Axi"'^^ , the 
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Fig. 4. The mean power spectra in the Parkes survey. The curves give the intensity 
Cj£ with and without source subtraction (full and dotted lines, respectively) and 
the polarization Cse and Cse (dash-dotted and dashed). 



logarithmic slopes vary with both the chosen submap and the range of £. 
In Table 1 we report angular-spectrum best-fit slopes for all of the submaps 
with no source subtraction. Columns 2 to 4 (primed quantities) refer to the 
range i < 500, and the remaining (unprimed quantities) to i < 800. Only for 
two submaps and for £ < 500 the intensity spectra are steeper than aj — 2; 
more often they are only moderately steeper than polarization spectra, or even 
hardly distinguishable from them (as far as slopes are concerned). 

We then have to combine the results from the 12 submaps in some way. A 
simple method is to compute weighted averages of the best-fit parameters; 
this implies weighting with statistical errors, with no regard to the strength 
of the emission in the submaps. Table 2 gives the mean parameters (including 
the normalization factors, denoted by Ax and A'x) obtained in this way, as 
well as the statistical errors; they are collected in the Table rows referring to 
averaging "method 1" . According to such results, the slope of the intensity 
spectrum slightly decreases for increasing i and can be hardly distinguished 
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Table 1 

Best-fit slopes of angular power spectra in the Parkes survey submaps with no source 
subtraction 



Submap 








Oil 


OlE 


OIB 


1 


1.79 


1.05 


1.24 


1.65 


1.28 


1.33 


2 


1.50 


1.28 


1.25 


1.64 


0.85 


1.12 


3 


2.54 


1.29 


1.45 


2.00 


1.31 


1.52 


4 


1.98 


1.74 


2.06 


1.46 


1.57 


1.74 


5 


2.87 


1.56 


1.47 


1.76 


1.49 


1.40 


6 


1.65 


1.16 


1.18 


1.89 


1.11 


0.96 


7 


1.56 


1.39 


1.78 


1.22 


1.56 


1.88 


8 


1.56 


2.26 


1.98 


1.43 


2.04 


1.90 


9 


0.74 


1.78 


2.02 


0.44 


1.56 


1.78 


10 


1.23 


1.28 


1.23 


1.02 


1.79 


1.25 


11 


1.97 


1.83 


1.61 


1.48 


1.48 


1.32 


12 


0.76 


2.09 


1.68 


0.96 


1.58 


1.45 



Table 2 

Best parameters for angular power spectra 
X Averaging A'-^ a'x Ax ax 





method 


Jy (beam area) ^ 




Jy (beam area) ^ 




I 


1 


241 ± 109 


1.71 ±0.18 


17 ±9 


1.37 ±0.13 


I 


2 


281?^ 


1.79 ±0.13 




1.6 ±0.13 


I 


2* 


{l.htf,^) ■ 10+5 


5.2 ±0.3 


{2.5tli) ■ 10+^ 


5.0 ±0.2 


I 


3 


12tf 


1.6 ±0.2 




1.36 ±0.13 


E 


1 


(8.5 ± 3.8) • 10-4 


1.57 ±0.11 


(11.5 ± 8.7) • 10-4 


1.44 ± 0.09 


E 


2 


{2.ltfl) . 10-4 


1.37 ±0.22 


(5.014;^) ■ 10-4 


1.53 ±0.11 


E 


3 


(1.31^:^) • 10-' 


1.31 ±0.18 


(2.7 ± 1.2) • 10-4 


1.45 ±0.1 


B 


1 


(16 ± 12) • 10-4 


1.58 ±0.10 


(4.4 ± 1.8) • 10-4 


1.46 ± 0.09 


B 


2 


(4.5li) • 10-4 


-, c;9+0.16 


(2.8l^i) • 10-' 




B 


3 


(8.2+^4) . 10-4 


1.35 ±0.20 


{Q.Stli) ■ 10-4 


1.32 ±0.13 



* After source subtraction 
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Fig. 5. The intensity and electric parity spectra compared to the fitting functions 
described in the text. The best-fit parameters refer to the range 36 < ^ < 800 and 
to averaging "method 2". 

from that of polarization spectra at large t. 



ai^aE^as^ 1.4 ^ 1.5, = 36 ^ 800, method 1). 



(5) 



Similar results are found with the simpler treatment with no signal cosine 
modulation; the differences with respect to the results in the Table are Actx — 
±0.1. 

It can be argued that because of the averaging method, the slopes given by 
Eq. (5) may be not representative of the full Parkes survey, since low-emission 
and high-emission subwindows arc considered on equal footing. An alternative 
approach is averaging the spectra of the 12 subwindows: Figure 4 reports the 
mean spectra, and includes magnetic parity and the intensity spectrum after 
source subtraction. We thereby performed best fits as above also on such 
mean spectra. This procedure, referred to as averaging "method 2" in Table 
2, clearly gives stronger weights to high-emission submaps. The results are 
quite different for the normalization factors (and especially for intensity), but 
no large variation is found for the slope factors. The slope of the intensity 
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spectrum however is slightly increased now, due to the scarce weight given to 
the nearly flat spectra of a couple of low-intensity submaps, the larger value 
being found in the smaller range £ = 36 -r- 500; no clear change however is 
found for polarization: 

q;j~1.8, q;e~q;s~1.4^1.5, = 36 ^ 500, method 2). (6) 

As shown by Fig. 5, the fitting functions with the parameters found in the 
larger range {£ < 800) describe very well the polarization spectra, whereas a 
modest discrepancy appears for the intensity spectrum at £ < 100. 

We also tried with a further procedure, renormahzing the spectra of the 12 
submaps to a same mean intensity value and then constructing mean curves. 
The parameters of such mean curves are referred to as averaging "method 3" 
in the Table, and show slightly smaller slopes. For the range £ = 36 -r- 800 

a/ ~ ~ as ~ 1.4, = 36 ^ 800, method 3). (7) 

A quite different result is found for the intensity maps with source subtraction. 
The slope is now much higher, aj ~ 5. At large £ the signal is very small, and 
it would be overcome by border effects in the absence of apodization even for 

e < 10^ 



3 Discussion 

The most attractive conclusion we can draw from the present work regards the 
angular spectrum slope of the polarized synchrotron emission up to ^ 10^. 
Since this slope substantially differs from that of polarized CMB, provided it 
can be extrapolated to the cosmological window the separation of cosmological 
and Galactic signals becomes easier. In particular, at a given frequency CMB 
polarization is expected to prevail better at smaller scales; this result is derived 
from Galactic-plane maps with no source subtraction. It is interesting to ob- 
serve that essentially the same slope, ~ ccb — 1-4, was found by Prunet et 
al. (1998) for polarized dust emission. We may speculate that this coincidence 
may be due to the similarity of polarizing/depolarizing mechanisms inscribed 
in the Galactic structure. 

On the other hand, our results on the total intensity spectrum might look 
somewhat intriguing at first sight. The original Parkes maps with no source 
subtraction support only a moderate slope, close to that of polarization spec- 
tra. This contrasts sharply with the results of Tegmark & Efstathiou (1996), 
Bouchet et al. (1996) and Bouchet & Gispert (1999); the discrepancy is much 
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Fig. 6. The angular spectra of polarized synchrotron extrapolated to 22, 32, and 60 

GHz (full line, from top to down) with the linear fits (dashed line), compared to 
the CMB E-spectrum computed for a standard CDM model (dot-dashed line) and 
a CDM with a secondary ionization optical depth Tion = 0.2 (dotted line). 

smaller with the analysis of Lasenby (1997) on the Tenerife patch. Removal 
of sources up to 2° wide, however, gives a much larger slope, beyond the more 
familiar result ctj 3. This fact can be readily interpreted, observing that 
in the Parkes survey the intensity angular spectrum at large i is dominated 
by Galactic-plane sources. An obvious consequence is that no real contradic- 
tion exists with previous work, since such results as a/ ~ 1.4 1.8 cannot be 
extrapolated outside the Galactic plane. 

We then should ask, which values of ax would be appropriate for the whole 
sky. Clearly the subtraction of 2° sources as described at the Parkes WEB site 
is a quite radical procedure, cutting out high-order harmonics very efficiently. 
It is then easy to infer that the full-sky aj must lie somewhere between 1.4 
and 5; however the existing support for the value a/ ~ 3 is not very strong, in 
view of the results of Lasenby (1997), supporting aj 2. On the other hand, 
we believe that our results on the polarization slope can be extrapolated out 
of the Galactic plane with some confidence. The inspection of Fig. 3 shows 
that for polarization no systematic differences of slope exist between high and 
low emission submaps. As already remarked, the overall normalization (rms 
signal) has much weaker variations than for intensity spectra. Thus the role 
of Galactic-plane sources is less important. Also, the agreement with Prunet 
et al. (1998) is impressive. 

Assuming that the polarization angular spectra are correct for the full sky, we 
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can try to extrapolate them to frequencies higher that 2.4 GHz and compare 
them to CMB spectra. Figure 6 shows such extrapolations to 22, 32 and 60 
GHz, where the average Parkes survey normalization and a synchrotron spec- 
tral index of 3 are assumed. Obviously this procedure is expected to overrate 
the synchrotron full-sky spectra. The comparison of synchrotron and CMB 
angular spectra shows that for £ > 10^ the latter prevail at 60 GHz in spite of 
the above Galactic-plane normalization. 

Finally, we may ask how this conclusion should be changed if the polarization 
slopes of Table 2 cannot apply to a full-sky survey. According to the reasoning 
expounded for intensity spectra, we can state that values ~ 1.4 can at least 
be used as lower limits to a^; and a^. Then the slope differences with respect 
to the CMB polarization spectra would still increase in comparison with the 
scenario described by Fig. 6, and the CMB signal would prevail even more 
strongly in the cosmological window. 
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